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Abstract 

We present an extension to a previous work to study the collapse of a radiating, slow-rotating self- 
gravitating relativistic configuration. In order to simulate dissipation effects due to the transfer of photons 
and/or neutrinos within the matter conhguration, we introduce the flux factor, the variable Eddington 
factor and a closure relation between them. Rotation in General Relativity is considered in the slow 
rotation approximation, i.e. tangential velocity of every fluid element is much less than the speed of light 
and the centrifugal forces are little compared with the gravitational ones. Solutions are properly matched, 
up to the first order in the Kerr parameter, to the exterior Kerr-Vaidya metric and the evolution of the 
physical variables are obtained inside the matter conhguration. To illustrate the method we explore 
the inhuence of the closure relations on the dynamics of three models with different equations of state 
and two functional form of the flux factor. We have found that, for the six closure relations considered, 
the matching conditions implies that a total diffusion regime can not be attained at the surface of the 
conhguration. It has also been obtained that the eccentricity at the surface of radiating conhgurations 
is greater for models near the diffusion approximation than for those in the free streaming out limit. At 
least for the static “seed” equations of state considered, the simulations we performed show that these 
models have differential rotation and that the more diffusive the model is, the slower it rotates. 
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1 Introduction 

Compact objects are one of the most fascinating objects known in our Universe. White dwarfs, neutron 
stars, quark stars, hyperon stars, hybrid stars and magnetars are thought to be relics from most of the cores 
of luminous stars which we believe to be born in supernova explosions. 

Core collapses are triggered by the implosion of the inner nucleus of a massive star {M^ ~ 8 — 2QMq) 
when its mass is in the limit of Chandrashekar {M^ore ~ 1.4Mq). During the implosion nearly all of an 
enormous gravitational binding energy /R ~ 5 x 10®^ ergs ~ 0.2Mc^) gained is stored as internal 

energy of a newly born, proto-neutron star (PNS) and driven by neutrino diffusion which cools this new type 
of compact object. Temporal and spectral characteristics of the neutrino emission depend on the rate at 
which they diffuse through the imploded PNS which, at this early stage, would have a mean density several 
times the standard nuclear density, p ~ 3M/ (47ri?^) « 7 x cm~^, with po — 2. x cm~^. The 

core density reaches up to (10 — 20)po during the cooling time tcooi ~ 5 to 10s, while the PNS de-leptonizes 
by neutrino emission, cools, contracts and spin-ups to form the final ultradense compact object [1, 2, 3, 4]. 

There is a consensus that the above standard scenario requires the description of General Relativity be¬ 
cause of the formidable gravitational fields arising during these processes. These powerful gravitational fields 
strongly couple hydrodynamics and neutrino flows within rotating matter configurations (see [5] for a good 
historical survey of previous works done at various levels on the problem of coupling the General Relativity, 
hydrodynamics, and radiation transport in spherical symmetry). Unfortunately despite a considerable effort 
that is been carried out by a significant group of people and institutions, presently we do not have any 
self-consistent model either analytical or numerical that includes all of those components in full details. 

Although an exterior metric of a rapidly rotating neutron agrees with the corresponding Kerr metric only 
to lowest order in the rotational velocity [6] , there have been many attempts to find a closed interior solution 
which matches smoothly to the gravitational field outside a rotating source (see [7, 8, 9] and references 
therein). In general these attempts have proved to be unsuccessful essentially because the considerable 
mathematical complexity in solving the Einstein equations [10]. It is only very recently, that there has been 
reported some progress in the analytical approach [11] which approximately match numerical solutions for 
rapidly rotating neutron stars [12, 13, 14] and has been used in studies of energy release [15, 16]. Recent 
numerical research has also considerably advanced our understanding of rotating relativistic stars [17]. There 
now exist several independent numerical codes for obtaining accurate models of rotating neutron stars in 
full General Relativity (see [7] and [8] for a good review on this subject). We can particularly mention a 
3D general-relativistic hydrodynamics code (GR Astro) [18, 19] and built from the Gactus Computational 
Toolkit [20] and some recent contribution for supernova scenario incorporating a better microphysics into 
rotating modeling of compact objects [21]. 

Simulations which include better microphysics in the form of realistic nuclear equations of state or 
neutrino transport have either been confined to spherical symmetry or restricted to newtonian gravity. 
Today, all available models analytical/numerical resembling some pieces of truth, demonstrate remarkable 
sensitivities to different physical aspects of the problem, in particular the treatment of neutrino transport 
and neutrino-matter interactions, the properties of the nuclear equation of state (EoS), multi-dimensional 
hydrodynamical processes, effects of rotation and general relativity. It is worth mentioning that is just 
recently, when it has become possible to obtain spherically symmetric general relativistic hydrodynamical 
core-collapse, treating the time and energy dependent neutrino transport in hydrodynamical simulations by 
considering a Boltzmann solver for the neutrino transport [22, 23], implementing multigroup flux-limited 
diffusion to Lagrangian Relativistic Hydrodynamics [5] or assuming the variable Eddington factor method 
to deal with the integro-differential character of the Boltzmann equation [24] . 

The present paper lies in between the traditional analytical and the emerging numerical descriptions of 
gravitational collapse. It follows a seminumerical approach which considers, under some general and reason¬ 
able physical assumptions, the evolution of a general relativistic rotating and radiating matter configurations. 
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The rationale behind this work is twofold, first it seems useful to consider relatively simple nonstatic mod¬ 
els to analyze some essential features of realistic situations that purely numerical solutions could hinder. 
Particularly, we will focus on the influence upon the evolution of matter configurations of the dissipation 
mechanism due to the emission of photons/neutrinos. Secondly it could be helpful for the evolving numerical 
codes to have testbed arena including General Relativity, rotation, dissipation and plausible EoS. 

The approach we follow to solve the Einstein Equations starts from heuristic assumptions relating density, 
pressure, radial matter velocity and choosing a known interior (analytical) static spherically symmetric ( 
considered as “seed”) solution to the Tolman-Oppenheimer-Volkov equation. This scheme transforms the 
Einstein partial differential equations into a system of ordinary differential equations for quantities evaluated 
at surfaces whose numerical solution, allows the modelling of the dynamics of the configuration. This method 
is an extension of the so called HJR [25], which has been successfully applied to a variety of astrophysical 
scenarios (see [26] and [27] and references therein) and which has been recently revisited [28, 29] in order to 
appreciate its intrinsic worth. 

We are going to consider the effects of rotation in General Relativity in the slow rotation approximation, 
i.e. up to the first order, thus we shall maintain only linear terms in the angular velocity of the local 
inertial frames. Thus, the effects of rotation are purely relativistic and manifest through the dragging of 
local inertial frames [30, 31]. This is understandable if we recall that in the newtonian theory where the 
parameter measuring the “strength” of rotation (the ratio of centrifugal acceleration to gravity at the equator) 
is not linear in the angular velocity but proportional to the square of it. The slow rotation approximation 
has recently proved to be very reliable for most astrophysical applications [13]. This assumption is very 
sensible because it considers that the tangential velocity of every fluid element is much less than the speed 
of light and the centrifugal forces are little compared with the gravitational ones. It is worth mentioning 
that the continuity of the first and the second fundamental {gij and Kij) forms across the matching surface 
are also fulfilled up to this order of approximation. 

Gonscious of the difficulties to cope with dissipation due to the emission of photons and/or neutrinos 
and, aware of the uncertainties of the microphysics when considering the interaction between radiation and 
ultradense matter, we extend a previous work [32] to study the collapse of a radiating, slow-rotating self- 
gravitating relativistic configuration by introducing a relation between the radiation energy flux density and 
the radiation energy density, i.e. the flux factor, f = J-/pn, and the so called variable Eddington factor, 
X = V/pr, relating the radiation pressure and the radiation energy density. We have also include a closure 
relation between both quantities, i.e., x = x(/)- the literature several closures have been introduced (see 
[33] for a comprehensive review and [34], [35] and [36] for more recent references) and most of them are 
consistent with the hiperbolicity and causality required by a relativistic theory [35]. 

With the above set of assumptions, i.e. seminumerical approach to solve the Einstein system, slow 
rotation approximation, and a particular closure relation between the flux and the variable Eddington factor, 
we shall explore the effect of dissipation on the evolution of the rotating radiating matter configuration. The 
outcomes from our simulations could represent rotating compact objects where the core rotates faster than 
the envelope. Therefore, the core can be supported by rapid rotation while the velocity of the fluid at the 
equator does not exceed the limit imposed by a fluid moving along a geodesic (the Kepler limit). Thus, 
differential rotation may play an important role for the stability of these remnants, since it can be very 
effective in increasing their maximum allowed mass. This effect was demonstrated in newtonian gravitation 
in [37] and was recently found by Shapiro and collaborators for general relativistic configurations having 
a polytropic EoS [38, 39]. For these rotating matter distribution we have found that boundary conditions 
imply that at the surface of the configuration, a total diffusion regime can not be attained. It has also been 
obtained that, with these coordinates, the eccentricity at the surface of radiating configurations, (up to first 
order) is greater for models near the diffusion limit approximation than for those in the free streaming out 
limit even though the rotation of configurations with dissipation near the diffusion limit appears to be slower 
than those near the free streaming out limit. We noticed that, at least for the static “seed” equations of state 
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considered, it seems that the lower the flux factor we have, the slower is the rotation of the configuration. 

The plan for the present work is the following. The next section contains an outline of the general con¬ 
ventions, notation used, the metric, the structure of the energy tensor and the corresponding held equations. 
Section 3 is devoted to describe the variable Eddington factor, the closure relations and the limits for the 
radiation held. Junction conditions and their consequences are considered in Section 4. The method is 
sketched in Section 5. We work out the modelling, previously studied for the spherical (nonrotating) case in 
Section 6. Finally some comments and conclusions are included in 7. 


2 Energy-momentum tensor and field equations 

2.1 The metric 


As in the previous work [32], let us consider a nonstatic, axially symmetric distribution of matter conformed 
by huid and radiation where the exterior metric, in radiation coordinates [40], is the Kerr-Vaidya metric[41]: 


ds^ =\l- 


2m (m) 


-I- o? cos^ ‘ 


du^ + 2dMdr — 2a sin^ Odrdcj) + 4a sin 9^ 


-(w) 


+ a^ cos^ 9 


dud(j) 


( 1 ) 


— -b cos‘‘ 9)d9^ — sin 


9 9 2m (u) r c? sin^ 9 

+OL ' ^ ' 


-b o? cos^ 9 


dcj) 


Here, m {u) is the total mass and a is the Kerr parameter representing angular momentum per unit mass 
in the weak held limit. It is worth mentioning at this point that the metric above is not a pure radiation 
solution and may be interpreted as such only asymptotically [42]. A pure rotating radiation solution may 
be found in reference [43]. However, as we shall show below although the interpretation of the Carmeli- 
Kaye metric in not completely clear, the model dependence of the considered effect is independent of the 
shape and the intensity of the emission pulse, and may be put in evidence even for a tiny radiated energy, 
AMrod = 10“^^M(0), which for any practical purpose corresponds to the Kerr metric[32]. The interior 
metric is written as [44] 


ds^ = dw^ -b 2d'iidr| — -b cos^ 9)d9‘^ -b 2de^^ sin^ 0 |l — — | dud(j> 

— 2e^^dsin^ 9drd(j) — sin^ 0 -b d^ -b 2d^ sin^ 


( 2 ) 


In the above equations (1) and (2), u = a:° is a time like coordinate, r = is the null coordinate and 
9 = and (j) = x^ are the usual angle coordinates. Local minkowskian coordinates {t,x,y,z) are related 
to Bondi radiation coordinates (it, r, 0, </>) by 


dt = ( \l—du + \l-^dr 


+ e^a sin^ ( 




(3) 


dx = w — (dr-b a sin^ 0d(()) ; dy = \/+ a cos^ 0d0 and dz = sin 0 -\/+ a cos^ 9d(j). (4) 


V 

The it-coordinate is the retarded time in flat space-time, therefore, it-constant surfaces are null cones open 
to the future. This last fact can be readily noticed from the relationships between the usual Schwarzschild 
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coordinates, (T, i?, 0,$), and Bondi’s radiation coordinates: 


u = T- 



0 = 0, r = R and 


<^= $ , 


(5) 


which are valid, at least, on the surface of the configuration. 

The Kerr parameter for the interior space-time (2) is denoted a and, for the present work it is relevant 
only (as well as a in eq. (1)) up to the first order. Notice that, in these coordinates, the r = rg = const, 
represent surfaces that are not spheres but oblate spheroids, whose eccentricity depends upon the interior 
Kerr parameter d and is given by 



( 6 ) 


Observe that this expression of eccentricity yielding the correct newtonian limit and corresponding to the 
natural definition in the context of metrics (1) and (2), is not invariantly defined . 

The metric elements fS and V in eq. (2), are functions of u, r and 9. A function m{u, r, 9) defined by 


V = e'^^ 


/ 2rh{u,r,9)r‘^ \ 

\ r^ + d? cos^ 9 ) 


(7) 


is the generalization, inside the distribution, of the “mass aspect” defined by Bondi and collaborators [45] 
and in the static limit coincides with the Schwarzschild mass. 


2.2 Energy-momentum tensor 

It is assumed that, for a local observer co-moving with a fluid having a velocity uj = (cu^, 0, oj^), the space-time 
contains: 

• an isotropic (pascalian) fluid represented by '' = diag {p, —P, —P, —P). Where p is the energy den¬ 
sity and P = Pr the radial pressure. Although the perfect pascalian fluid assumption (i.e. Pr = P±) is 
supported by solid observational and theoretical grounds, an increasing amount of theoretical evidence 
strongly suggests that, for certain density ranges, a variety of very interesting physical phenomena may 
take place giving rise to local anisotropy (see [46] and references therein). 

• a radiation field of specific intensity I(r’, t; n, v) given through 

df = 1(7’, t; n, z/)dS' cos (/j dO du dt, (8) 

with ip the angle between n and the normal to dS" and where d£ is defined as the energy transported 
by a radiation of frequencies (y, v y dw) in time dt, crossing a surface element dS, through the solid 
angle around ft, i.e. dO = sin0d0d7/7 = —d^dif. 

As in classical radiative transfer theory, for a planar geometry the moments of I(r, t; n, v) can be written 
as [47, 48, 49] 

^ nOO f~^ 1 

Pfl = 2 y J dp l{r,t,fi,n), ^ 2 J J 9 -^'r,t,fi,v) (9) 

and ^ 

7^ = - / dn j dp p^Hr,t-,n,v) . (10) 

2 Jq 

Physically, pp> , T and V, represent the radiation contribution to the: energy density, energy flux 
density and radial pressure, respectively. 
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From the above assumptions the energy momentum tensor can be written as where the 

material part is and the corresponding term for the radiation field, can be written as [47, 48]: 


-TO 0 \ 

VO 0 

0 Upr-V) 0 

0 0 \{PR-V) ) 

Notice the induced anisotropy in the due to the radiation field. 

Then, the energy-momentum tensor in the local co-moving frame takes the following form: 


rpR _ 

^ fliy 


/ PR 
-T 
0 

V 0 


( 11 ) 


= 


ip + Pr) + P + 2 iPR ~ 


~ + o iP^ ~ ^) ) o(^^ ~ PR)XtiX'' + Pfj.Ui' + FvUfi 


where = diag [I,-!,-!,-!) = (1,0, 0,0), = (0,1,0,0) and = (0,-J^,0,0). 

Now following [32] , in order to find the energy momentum tensor as seen by this observer co-moving with 
the fluid, we should perform an infinitesimal rotation around the symmetry axis, i.e. 


( P + PR 
-T 
0 

\ \V{ipR-V) 


-T 0 

P + V 0 

0 P + \{pR-V) 
-VT 0 


lV{3pR-V) \ 
-VT 


0 

P+\iPR-V) / 


( 12 ) 


where V {u, r, 9) is associated with the local “dragging of inertial frames” effect, which in the slow rotation 
limit V will also be taken up to first order. 

Once minkowskian co-moving energy momentum tensor is built in terms of physical observables on a local 
frame (p, P, pR , T, 7^,and I?), it can be transformed from the local minkowskian co-moving coordinates 
(t, a;, y, z) to the curvilinear not co-moving Bondi coordinates (u, r, 9, (p) as 


where L\{u!i) is a Lorentz boost, written as 


(13) 


( ^ 




-lUJx 

0 


1-b 


T^x 


V -7^2 


-- 


-l^z 
LJx^z ( 7 ~ 1 ) 


1-b 


0 


with 7 = 




and 


2 2,2 




(14) 


Observe that d3P/dx°‘ are coordinate transformations connecting {t, x, y, z) with (tt, r, 9, 4>) which can be 
identify from equations (3) through (4). 

In radiation coordinates the radial and orbital velocities of matter are given by 


dr V ujx 

du r 1 — ^ 2 : 


and 


dcp 

du 


_ 1 ^0 

1 — LOx r sin (0) V r ’ 


(15) 


respectively. 

Now, using the metric (2), the energy momentum tensor (12), the transformation (13) and considering 
the slow rotation limit (i.e. first order in the orbital velocity uoz, Kerr parameter d and the dragging 
function 27), we can write the Einstein (see Appendix 9.1). As in reference [32] only six of the eight physical 
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variables P, P, Pr,^ and V), can be algebraically obtained, in terms of the metric functions 

(3{u,r,9) and m{u,r^9) and their derivatives, from field equations (69) through (77) . Therefore, more 
information (equations) has to be provided to this system in order to solve the physical variables. The idea 
will be to supply relations among the radiation physical variables PR,P, and V. Next section will be devoted 
to describe these essential relations. 


3 Closures relations and the limits for the radiation field 

In order to deal with more realistic scenarios, the microphysical framework of the interrelation between matter 
and radiation have to be considered. The relativistic Boltzmann Transport Equation must be coupled to 
the hydrodynamic equations in order to obtain the evolution of the system as well as the spectrum and 
angular distribution of the radiation field [47] Neglecting effects as polarization, dispersion and coherence, a 
covariant special relativistic equation of radiation transport has been proposed as [50, 51]: 

(u^ + n i + u, + 

where = 1 with = 0; the four velocity of the fluid is p is proper density of the medium; the 
quantities eg and k are the emissivity and the absorption coefficient, respectively. This transfer equation has 
several important difficulties. The most important are: the lack of information about the coupling between 
radiation and ultradense matter and its mathematical complexity, although some understanding is emerging 
recently [52, 53]. 

One of the possible strategies to circumvent the difficulty of solving the radiation transfer equation is 
to consider one of the two physical reasonable limits for the radiation field which describe a significant 
variety of astrophysical scenarios [48]. The free streaming out limit assumes that radiation (neutrinos and/or 
photons) mean free path is of the order of the dimension of the sphere. This was the case considered in [32] 
and it can be expressed as 

PR=T = V = e. (17) 

The other limit for the radiation field is the diffusion limit approximation, where radiation is considered to 
flow with a mean free path much smaller than the characteristic length of the system. Within this limit, 
radiation is locally isotropic and we have 

Pr = iV and T = q . (18) 

In order to simulate more realistically the matter and radiation interaction, it seems more reasonable 
to have a parameter which varies between the above mentioned limits. This is the idea of the flux and the 
variable Eddington factor an they can be summarized as follows. From equations (9) through (10) it is 
convenient to define the following normalized quantities 

(p{f,t,fl) = f = / p{r,t,n)ndfl and K = / (p{r,t,VL)fi®ndVt, (19) 

Pr a Att J Att 

The Eddington factor is, defined as the eigenvalue of the Pressure Tensor corresponding to the eigenvector 
n (unitary vector in the direction of the energy flux ), i.e. Kjn^ = XU* [54]. Thus, 

f ^ = /n* and K i {(1 - + (3% - } • (20) 

In the one-dimensional case the above equations lead to 

T V 

f = - and X = — ■ (21) 

PR PR 
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Which are called the flux and the variable Eddington factor, respectively. 

In order to “close” this problem and to algebraically obtain six of the above mentioned physical variables, 
namely Wz, p, P, T, and I?, from field equations (69) through (77) and the radiation parameter (21) (or 
in general (19)) we need to state a relation between /, and y. It is easy to perceive that such a relation 
could exist. In fact, it is noticeable that in the corresponding limits for the radiation field, i.e. dijfusion limit 
approximation and free streaming out we have 

P = \PR ^ f —>0andx=i'| . 

;^0</<l and -< xif) < 1 (22) 

P=V =pR ^ / = 1 and X = 1 J 

Causality requirement implies the following supplementary conditions on / and x^ order to define a 
physically plausible region in the {/, Xj^x/d/} space [35] 

||/||<1, f<X<l and (23) 


Closure 

xif) 

dx 

d/ f=i 

dy 

.f=0 

Lorentz-Eddington (LE) 


2 

0 

Bowers- Wilson 

|(l-/ + 3r) 

5 

3 

1 

3 

Janka (Monte Carlo) (MC) 


2.28 

0 

Maximum Packing (MP) 

4(l-2/ + 4>) 

2 

'2 

3 

Minerbo (Mi) 

X (/) = 1 — 2^ where / = coth k — ^ 

2 

0 

Levermore-Pomraning 

xif) = f coth /3 where / = coth /3 - | 

1 

0 


Table 1: Closure Relations and some of their physical acceptability conditions 

There are several closure relations reported in the literature (see two recent comprehensive discussions on 
this subject in [35] and [36] and references therein). Few of them are simply ad hoc relations that smoothly 
interpolate the radiation field between the diffusive and free-streaming regimes. Others, are derived from a 
maximum entropy principle or from a given, or assumed, angular dependence of the radiative distribution 
functions. Even one of them has been motivated from direct transport calculations. Six of the most frequent 
found closure relations are listed in Table 1. In this list the first four could be considered as “analytical” 
closure relations, while the last two are referred as numerical, because for a given flux factor /, the nonlinear 
equation / = coth/3— (1//3) has to be numerically solved in order to obtain the variable Eddington factor x- 

We are going to explore some of the effects of dissipation on the evolution of slowing rotating radiating 
matter configuration in General Relativity. We shall evaluate how independent are these effects from an 
explicit closure relation and/or a specific EoS chosen. Particularly, some results concerning the influence 
of the junction conditions on the eccentricity and the radiation scheme evaluated at the surface, will be 
presented in the next section. The strategy we follow to close the system of Einstein field equations with a 
radiation field, contrasts with the standard iterative method for solving the moment equations (9), and (10), 
starting from an estimated Eddington factor (see [48] and [24] and references therein). 

Thus, by using equations (21), the Einstein field equations (69)-(77) can be re-written in terms of the flux 
and the Eddington factors (see Appendix 9.2) . Again, in principle, for all cases listed in Table 1 (including 
numerical closures relations) it is possible to obtain the remaining six physical variables, uJx,tUz, P, P,P, and 
V, from the system (78) through (86), in terms of the Kerr parameter, a, the flux factor, / , the metric 
functions l3{u,r,0),rh{u,r,6) and their derivatives. 
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4 Eddington factor, flux factor and junction conditions 

In this section, following [32] we should match the interior fluid spheroid to the exterior Kerr-Vaidya solution 
(equation (1)). Therefore, the continuity of the first and the second fundamental [gij and Kij) forms across 
the matching surface are needed. These requirements are equivalent to demand the continuity of the tetrad 
components and spin coefficient of the metrics (1) and (2) across the boundary surface r = a{u) [55]. 


4.1 Junction conditions for a slowly rotating configuration 

The Newman-Penrose null tetrad components for the metrics (1) and (2) (see reference [32] to find the 
expressions of the spin coefficients for these two metrics) are for the exterior metric 


It^ = Sft; n^=Si^- 


1 


1 - 


2mr 


r^ + o? cos^ iff) 




ia sin 6 {S!^ — S!^) + Sg + i esc 0(5^ 


and 


11^ = n>^ = S^- -e-^^ 

r I “2 


1 - 


2rhr 


r^ + OL^ cos^ 0 




V2{r + ia cos 9) 

ia sin 9 {S!^ — S^) + Sg +i esc 96^ 


V2{r + ia cos 9) 


(24) 


(25) 


for the interior metric ( 2 ) 

The continuity of the tetrad components across the boundary surface r = a(u) implies 

Pa = 0 ; rha = m and 5a = a , 
and the continuity of the spin coefficients r, 7 , and v lead to 

/ 2777- \ 

/3la 1 -- /3oa = Pla = m 2 a = 0 and a {Pia - Poa) = a (moa - rriQa + rhia) = 0, 

\ a J 2a 


(26) 


(27) 


which means that {Pia — Poa) and (moo — rhoa + rhia) are of order a. 

Now, evaluating the field equations (69) - (77) (or equivalently (78) through ( 86 ), ) at r = a{u) and 
considering the above results (26), we obtain that, on these coordinates and up to the first order in a, the 
metric coefficients P and rh are independent of the angular variable and consequently the physical variables: 
uJx, p, P, and T are also d—independent (see [32] for details). 

Next, expanding P near the surface, Poa + a/3ia = 0, in the two first equations (27) and using that P is 
continuous and vanishes at the outside the matter configuration we obtain that 




d = I 1 — 


2.fha 


{Pa + PRa ~ Pa) ^xa ~ {Pa + Pa ~ Pa) 


_ {Pa + PRa ~\~ Pa~\~ Pa ~ a) (1 ~ ^xa) _ 


where a = dr/du. On the other hand, from equation (15) it follows that 


d = 1 - 


2fha 




1 UJx 


(28) 


(29) 


Equating (28) and (29), it is obtained that the emerging energy flux density compensates the total pressure 
(hydrodynamic and radiation) inside the configuration [56], i.e. 


Pa — Pa P Pa ■ 


(30) 
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Now, expanding (3 and m near the surface in the third equation (27), we conclude that 


/3ia (1 + a) « a Pia (1 + a) = v{u)a and mia (1 + a) « a mia (1 + a) = q{u)a . (31) 


Where v{u) and q{u) are arbitrary functions of the time-like coordinate u with |u(u)| ^ 1 and \q{u)\ ^ 1 in 
order to keep valid the approximation. Additionally, by using field equations (70) and (71), we get that 


/3ia (1 -b a) « a 
and 

rhia (1 -b a) « a 


27ra I 1 — uj-r 


2 rha 


2 TJla 1 /. \ 

1-I (1 + ^xa) 


{Pa ~b PRa ^ a] ~ 


(32) 


a^dTT 
1 + 


1 ^xa ~ I \Pa “b PRa -F a] ~ O , 


which impose restrictions on the physical variable evaluated at the surface of the distribution. 
Also, it follows from the junction conditions (27) that 


l-ba-^^^=TOia 2a/3ia ^1 - = mia (1 - Wa;a) • (33) 

Thus, we obtain an expression relating v{u) and <?(«), namely 

2a[l + d- v(^u) = q{u) —r fl- v{u) = q{u). (34) 

\ a J (1 - LOooa) V a / 

However, it can be checked by simple inspection that, because neither the held equations nor the junction 
conditions impose further limitations on these functions of u, one of them remains completely arbitrary for 
each model. 

More over, at least for some models we have 



2a 

(1 ^xa) 



1 , 


(35) 


which becomes useful when selecting the order of magnitude of the initial parameters (rha, a and ujxa) for 
the modelling of slowly rotating collapsing conhgurations (i.e. hrst order in the orbital velocity tOz, Kerr 
parameter a and the dragging function V), worked out in section 6.3. 

The next section will be devoted to explore consequences of the radiation held, T^, in collapsing conhg¬ 
urations in the slow rotation approximation. 


4.2 The limits for the flux factor and the eccentricity 

Now, from equations (21), and (30), we have 

Pa = (l - ^) ^Pa= Fa ' (36) 

If we assume that the hydrodynamic pressure and the outgoing energy hux have to be positive then we have 
that (1 — (xalfa)) > 0. Figure 1 displays this factor for the diherent closure relations in Table 1. It is clear 
from this hgure that the boundary conditions compel the impossibility to attain total diffusion regime ( 
i.e. / = 0) at the surface of the conhguration. The roots of each curve representing a closure relation dehne 
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Closure 

fr=a 

e 

Lorentz-Eddington 

f < fLE\r=a ^ ^ 

A<«*<a(i + 1 :L( 

Bowers- Wilson 

5 < fBw\r=a < 1 

A < cbw < a ^1 -|- 2— 

Janka (Monte Carlo) 

0.39 < fMcl=a < 1 

A < CMC < A (^1-h 1.545^ j 

Maximum Packing 

J < fMp\r=a ^ 1 

A < cmp ^ a ^1 -|- 3—^ 

Minerbo 

0.40 < fM\r=a < 1 

A < CM < A 1^1 -h 1.488^ j 

Levermore-Pomraning 

fLp\r^a = 1 

CLP = A 


Table 2: Limits for the flux factors and the eccentricity for the different closure relations 


the interval of acceptability for the values of the flux factor / . These intervals are displayed in the second 
column of Table 2. As it can be appreciated form Figure 1, only the Levermore-Pomraning closure relation 
does not meet this requirement. Up to the precision of our numerical calculation, the acceptable value 
for / that guarantees the positiveness of the hydrodynamic pressure is 1. Thus, considering the Levermore- 
Pomraning closure relation, junction conditions, up to the first order in Kerr rotation parameter, only allows 
free streaming out at the surface for this slowly rotating matter distribution. On the other hand, Minerho 
closure relation seems to admits transport mechanism closer to the diffusion limit. 

Because all these results emerge from the junction conditions that couple the internal and the external 
solutions, they are valid not only for axisymmetric configurations but also for spherical ones. It is also 
independent of the EoS and is present for all the closure relations we have listed in the Table 1. 

Finally, expanding ( 6 ) for d ^ 1 , we get d[57] 


1 - 11 


e = —a- -—a 
r, 2 rj 


(37) 


as expected, up to first order, the eccentricity is proportional to a. Now, using the held equation (70) 
evaluated at the surface r = a{u), (31) and (34) we are lead to 


On, = o = 


27ra [ Pa + Pa 


(1 - fa) 

fa 


1 - 


2rha 


(1 + i^xa) ( 1 ~ 


2mc 


v{u) 


(38) 


and the surface eccentricity can be re-written as 


27r 


— 


Pa P Pc 


(1 - fa) 


fa 


1 - 


2rha 


(1 -|- UJxa) ( 1 ~ 


2fnr. 


v(u) 


(39) 


Notice that, v{u) remains completely arbitrary and its choice completes the characterization of the model. 
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Because of the range of acceptability for the flux factor we also obtain a range for the eccentricity, i.e. 

2ma 


/minx</<l ^A<e„<A 1 + 


(1 - fa) 
fa 


where A = 


27rpa 1 - W; 


(1 + UJxa) ( 1 ~ 


2fria 


(40) 


v{u) 


Therefore, it is clear from Table 2 that radiation mechanism affects the oblateness of the configuration. This 
is to say, in these coordinates up to first order in da and for those models having A > 0, the eccentricity at 
the surface of a radiating configuration is greater for models near the diffusion limit approximation than for 
those in the free streaming out limit. Again, this result is also valid for any EoS with 0 ^ v{u) ^ 1 and it is 
present for all the closure relations in the Table 1. 

In the next section we shall explore the effect of these closure relations on the collapse of a radiating, 
slow-rotating self-gravitating relativistic configuration. 


5 The HJR method and the surface equations 


In order to obtain the evolution of the profiles of the physical variables, Ux, P, P, P, and V, we use an 
extension of the HJR method [25] to axially symmetric slowly rotating case [32]. 

First, we define two auxiliary variables which, in terms of the Eddington and the Flux factor, can be 
written as 

p -f pr — J-—uJx{P -\-P — _ P ~ ^xP + j(l ~ ~ ^x{x ~ f ))P 


p = _ 

1 + 1 + ^x 

and is called the effective density and, correspondingly, the effective pressure is 


P = 


P + V- P-OJx {p + PR-P) _ P - ^xP+ jix - f - ^xi^ - f))P 


(41) 


(42) 


1 + Ulx 1 + iOx 

With these effective variables, the metric elements (equations (70) and (71)) can be formally integrated as 


B{u,r)= / 2TTr-, - rrz^dT and rhiu.r) = / 4tt r'^ pdf . 

Ja (1-^) Jo 


(43) 


Thus, if the r dependence of P and p are known, we can get the metric functions fh and /3 up to some 
functions of u related to the boundary conditions. This is one of the key points to transform the Einstein 
System into a system of (coupled nonlinear) ordinary differential equations on the time-like coordinate. 
Physically, the rationale behind the assumption on the r dependence of the effective variables P and p, 
can be grasped in terms of the characteristic times for different processes involved in a collapse scenario. If 
the hydrostatic time scale Thy dr, which is of the order ~ Xj \JGp (where G is the gravitational constant 
and p denotes the mean density) is much smaller than the Kelvin-Helmholtz time scale ( Tkh ), then in a 
first approximation the inertial terms in the equation of motion can be ignored [3]. Therefore in this first 
approximation (quasi-stationary approximation) the r dependence of P and p are the same as in the static 
solution. Then the assumption that the effective variables (41) and (42) have the same r dependence as the 
physical variables of the static situation, represents a correction to that approximation, and is expected to 
yield good results whenever Tkh ^ Thy dr- Fortunately enough, Tkh ^ Thy dr, for almost all kind of 
stellar objects. Recently this rationale becomes intelligible and finds full justification within the context of 
a suitable definition of the post-quasi-static approximation for the gravitational collapse [28, 29]. 
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Those functions of the time-like coordinate u that remain arbitrary can be obtained from a system of 
ordinary differential equations {The System of Surface Equations, SSE) emerging from the junction condi¬ 
tions and both the field equations and some kinematic definitions evaluated at the boundary surface. The 
first surface equation is (29): 

A = F{n- 1 ) . (44) 

Where we have scaled the radius a, the total mass rha = m and the timelike coordinate u by the total 
initial mass, m{u = 0) = m(0), i.e. 


A = 


a 

m{0) ’ 


M = 


m(0) ’ 


u = 


to ( 0 ) 


(45) 


and we have also defined 

F = 1 — and fl = -— - . (46) 

Again, the dot over the variable represents the derivative with respect to the time-like coordinate. The 
second Surface Equation emerges from the evaluation of equation (78) at r = o+q. It takes the form of 


M = -FL , 


(47) 


where L representing the total luminosity can be written as 

L = AirA^Ta (20 - 1). (48) 

Now, using above equation (44) and definitions (45) and (46); we can re-state equation (47) as 


F 

F 


2L + {1 - F){n - 1) 


(49) 


Finally, after some straightforward manipulations, starting from held equations (79), (80) and (81), it is 
obtained 


. 2/3 


p + P 

_ 2m 


^ p + p 
dr 


1 - 


AirrP - 


-2 

r 


P+\{PR-V)-P], 


(50) 


which is the generalization of Tolman-Oppenheimer-Volkov (TOV) equation for any dynamic radiative sit¬ 
uation. The third Surface Equation can be obtained evaluating (50) at r = a+o, and it takes the form of: 


0 = ^ \ ^ \ I 

Fl F Pa Pa Apa 


P. + ^^(x (/)-!) )+(0-l) 


F^Pla 

Pa 


where 


R = 


dP 

dr 


p + P 

_ 2m 


AnrP ' 


m 

72 


47rA(l-30)pa S + F\ 

O J 

(51) 

(52) 


Equations (44), (49) and (51) conform the SSE which coincides with the spherically symmetric case [56] 
because that up to the Hrst order in d, the metric functions m and /3 are found to be independent on the 
angular variables. This system may be integrated numerically for any given radial dependence of the effective 
variables, providing the total luminosity, a closure relation and a hux factor /. The remaining two equations 
(85) and (86) provide a simple 0-dependence on the physical variables Uz, and V, i.e. 


ujz = d sin 9 y [rh, /?; their derivatives; r] and V = d sin 9 Z [rh, /3; their derivatives; r] (53) 
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The restriction due to the junctions conditions for slowly rotating spheroids (32) can be re-written in 
terms of the effective variables as 

2'kA ^°' ^ -I- = V (u) a and 47r 'pa + -^ = q {u) a, (54) 

and equation, at least for some models, (35) can be re-phrased in a very compact form: 

2AFn ~ 1 (55) 

Again, this equation becomes very useful when selecting a set of initial conditions to integrate the SSE. 

For completeness, we outline here a brief resume of the HJR method for isotropic slowly rotating radiating 
fluid spheres (see [32], for details): 

1. Take a static interior solution of the Einstein Equations for a fluid with spherical symmetry, pstatic = 
p{r) and PstaUc = P{r). 

2. Assume that the r dependence of P and p are the same as that of PstaUc and pstatic, respectively. Be 
aware of the boundary condition: 

Pa = -^xaPa- (56) 

and equations (54). 

3. With the r dependence of P and p and using (43), we get metric elements rh and [3 up to some functions 
of u. 

4. In order to obtain these unknown functions of u, we integrate SSE: (44), (49) and (51). The first two, 
equations (44) and (49), are model independent, and the third one, (51), depends of the particular 
choice of the EoS. 

5. One has four unknown functions of u for the SSE. These functions are: boundary radius A, the velocity 
of the boundary surface (related to O), the total mass M (related to F) and the “total luminosity” L 
Providing one of these functions, a closure relation and the flux factor /, the SSE can be integrated 
for any particular set of initial data a that fulfill equation (55). 

6 . By substituting the result of the integration in the expressions for ifi and /3, these metric functions 
become completely determined. 

7. Again, once we have provided a closure relation and the flux factor /, the set of matter variables, 
ujx, P, P, and T can be algebraically found for any part of the sphere by using the field equations 
(78)-(81); rotational physical variables, lOz and I?, can be obtained from the remaining significant, two 
held equations (85) and (86). Finally, radiations variables, pn and V, emerge from (21), introducing 
the fiux factor, /, the variable Eddington factor y and any closure relation. 

6 Modelling slowly rotating matter configurations 

In order to explore the influence of the dissipation mechanism and the effect of closure relation on the 
gravitational collapse of slowly rotating matter configurations, we shall work out three models previously 
studied for spherical (nonrotating) cases. We shall work out three different EoS: Schwarzschild-like [25, 58], 
Tolman IV-like [56, 58, 59] and Tolman VFlike [25, 56, 58]. 
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6.1 The models 


The first family of solutions to be considered is the slowly rotating Schwarzschild-like model. In the static 
limit this model represents an incompressible fluid with constant density. It is the same example presented 
in ref. [32] but for the present case we have included the flux & Eddington factors (21) and a closure relation 
from Table 1. The corresponding effective density and pressure can be written as 


p = k[u) = 


3 1-F 
Stt 


and P = k{u) 


ig{u) [l — {8Tr/3)k{u)P 


, 1/2 


- 1 


3 — 8g(u) [1 — (87r/3)fc(u)r^] 


211/2 


where the function g{u) can be determined from the boundary condition (56), as 

3-20 


giu) = 


[1 — (87r/3)fc(u)a^]^^^ 


Third surface equation (51) is 


n = 


-n 

1 - F 


3(1-F)2(2n-l)(fl-1) F 


2An 


F 


(57) 


(58) 


(59) 


The second EoS to be discussed corresponds to the slowly rotating Tolman-IV-like model. This model 
exhibits, in the static limit at the center, the EoS for pure radiation, i.e. P/p ~ 1/3. The effective density 
and pressure for this case can be expressed as 


1 + 3 


P = 


Z{u) 

W(u) 


’Wiu) 


1 - 


W{u) 


8TrZ{u) 


1 + 2 


-Zlu) 


where 

Z{u) = - 
with 


A^ [7(1 -F)+ 2n(F -2)-p] 


(l + 2^) ^ 


and W{u) = 


and P = 


1 _ Z(n) _ o p 
W{u) '^W(u) 


1 + 2 


(60) 




(F(i + 20)-l + 77) 


2(F-l)(2fl + 3) ’ 2 [E(2-3n)++2(60-5) + ! +(F-1)?7]^ 

p = Vl + F{22 - 20n) + +2(402 + 20 - 23). 


equation 

For this model, the third surface equation can be written as 

O = 0i + $+ + r 


(61) 


where the expression for the coefficients 0,5*, and T in terms of the Surface Variables and their derivatives 
(i.e. A, F, O, A,and +) are sketched in the Appendix 9.3. 

The third family of models is inspired on Tolman VI static solution, which approaches the one of a highly 
relativistic Fermi gas, with the corresponding adiabatic exponent of 4/3. For this case we have 


P = 


3h(u) 


1 


87rr2 


( 1 -+) 


and 


+ = 


h{u) (1 — 9d{u)r) 
r2 (1 — d{u)r) 


As before, the function d{u) is determined from the equation (56), thus, we obtain 

M \ 1 411-1 

d{u) = -- 


3 ^(4^ - 3) ■ 


(62) 


(63) 
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We find that the third surface equation is 


4F2(F - 1)^3(2172 _ + 2) + 2FAn{n - 1)^ + F(3 - 2F - F^)n^ - F{F - 1)2(1 - 30) 

2AF(F- l)(n- 1)2 


(64) 


6.2 General considerations 


For all the models, we have selected a set of initial conditions and physical parameters that resemble, as 
much as possible, interesting astrophysical scenarios. We have chosen 


• Typical values of these conditions that resembles young neutron stars. Notice that because the coupling 
restriction for slow rotation assumption (35) (or its equivalent in the adimensional variables (55)) the 
initial radius of the configuration becomes 4 times greater than the typical neutron star radius. 

• In our simulations, we have imposed that the energy conditions for perfect fluid be satisfied. In addition, 

the restrictions —1 < < 1 and r > 2rh{u,r) at any shell within the matter configuration are 

also fulfilled. 


• As it was pointed out at the end of the preceding section, the evolution of one of the variable at the 
surface (boundary radius A, the velocity of the boundary surface (related to O), the total mass M 
(related to F) or the “total luminosity” L) has to be provided. For the present simulation the evolution 
of the luminosity profile, T(u), is given as a Gaussian pulse centered a,t u = Up 


-M = L = 


AM, 


rad 


1 f u — Ur 




where A is the width of the pulse and AM^ad is the total mass lost in the process. 


(65) 


• Throughout the simulations equations (31) (or equivalently (35) or (55)) are constantly checked in 
order to verify the validity of the approximation. 


6.3 Modeling radiation transfer scenarios 

We would like to explore how dissipation affects the dynamics of these three types of slowly rotating matter 
distributions. For each of the above EoS and the several closure relations listed in Table 1, we shall work 
out simulations with: 

1. Matter configurations with constant flux factor f. We study several radiation transfer environ¬ 
ments ranging from the collapse of opaque matter distribution where / = 0.426 (close to a diffusion 
regime) to more transparent matter configuration where / = 0.930 where the radiation transport 
mechanism is described near the free streaming out limit approximation. 

2. Matter configuration with variable flux factor f = f{r). For this case we study the effect of a 
variable flux profile as 


/ 



p-Cioot-x) f < f f 

^ Jcore \ Jsurja 

1 -I- 


( 66 ) 


on the orbital velocity at the equator. We have defined fcore as the flux factor at the inner core and 
fsurf ace the flux factor at the surface of the distribution. The parameters Xt = rt/m{0) represents 
the cutoff region where the transition of the dissipation mechanism takes place and C regulates how 
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sharp or smooth is the transition between the flux factors at two regions (see Figure 2). The idea with 
this variable flux factor is to allow conflgurations with more opaque matter at the inner core and more 
transparent mass shells at the outer mantle. In many Astrophysical scenarios, radiation diffuses out 
from a central opaque region (x(/ = 0) = 1/3) to the transparent boundary (x(/ = 1) = 1). 

6.4 Parameters and astrophysical scenarios 

Concerning the luminosity (equation (65)), all our models have been been simulated with 

/^Mrad = 2.00 X 10"“ M(0), A = 0.74 x W^s, tp = 1.48 x lO-^s. (67) 

For the modelling with the Lorentz-Eddington closure relation and the three above mentioned EoS we 
have set initial conditions and parameters to have the following values: 


m(0) = 

I.OMq 

a = 

10"^ 

4(0) = 

6,000 

d (0) = T'gurface — 

44, 500 m 

U(0) = 

0.999 

(0) = 

-0.00101c 


The Lorentz-Eddington closure relation was initially proposed by C.D. Levermore in the early 80’s [33] on 
the basis of geometrical considerations, for the case of stationary medium. Later, it has been reobtained 
by other authors from different perspectives, i.e., thermodynamical point of view with maximum entropy 
principles [54] and Information Theory with the energy flux taken as a constraint [34] . 

For all other closure relations listed in Table 1, the set of initial conditions and parameters for the 
modelling performed with the above seed EoS are: 


m(0) = 

I.OMq 

a = 10“^ 

4(0) = 

41,104 

^ a (0) = r Surf ace = 36, 500 m 

U(0) = 

0.999 

^a;a(0)= -0.00101c 


Concerning the flux factor, /, we have considered two different scenarios for the three equations of state 
and for the six closure relations. First, for simplicity, we are going to consider the case / = const. The second, 
and more realistic, scenario, described by equation (66) and displayed in Figure 2, has the corresponding the 
set of parameters: 


/surface = ^; /core = 0.902; ( = 10 and Tj = 14,833m (68) 

7 Summary of results, comments and conclusions 

We have extended a previous work [32] to study the collapse of a radiating, slowly rotating self-gravitating 
relativistic configuration by introducing the flux factor., the Variable Eddington Factor (equations (21)) and 
a closure relation (Table 1). With this extension, in principle, it is possible to implement the seminumerical 
approach to simulate a variety of radiation transport mechanism within a slowly rotating radiating matter 
distribution that can be used as evaluation testbeds for emerging full-numerical environments. 

Using this seminumerical approach it has been explored the influence of the closure relation on the 
dynamics of collapsing slowly rotating relativistic bodies. The results we found can be summarized as 
follows. First, it has been obtained that, for all closure relations listed in Table 1, the junction conditions 
imply that total diffusion regime can not be attained at the surface of the configuration. This result is exact 
and general, independent of the EoS and valid for spherical and axisymmetric matter configurations. 

Secondly (also related to the junction conditions), it is evident from equation (39), that the eccentricity 
at the surface of radiating conflgurations (up to first order in a) is greater for models near the diffusion limit 


Aguirre, Nunez and Soldovieri 


17 












Variable Eddington Factor and Radiating Slowly Rotating... 


approximation than for those in the free streaming out limit. Again, this result is EoS independent and is 
present for all closure relations we have studied. 

The third result is the significant influence of the closure relations on the dynamic of physical variables. 
For the more realistic scenario concerning a variable flux factor (equation (66)) and for the first order 
approximation in rotating parameters, we have found qualitatively different behavior for the profiles of the 
tangential velocity As it can be appreciated from figures 3-6 we obtain models having outer layers 
rotating faster than the inner ones or experimenting counter rotation for some mass shells surrounding the 
nucleus. 

Finally some other results emerge from our modeling. Concerning newtonian theory, it is clear that the 
gravitational effects of rotation are purely relativistic and as it can be understood recalling that the newto¬ 
nian parameter measuring the “strength” of rotation is not linear in the angular velocity but proportional 
to the square of it. In this work, rotation is considered in the slow approximation limit, i.e., we dealing 
with situations where the tangential velocity of every fluid element is much less than the speed of light and 
the centrifugal forces are little compared with the gravitational ones. Despite that the physical rationale 
to study the gravitational effects of rotation in the slow rotation approximation seems to be justified, this 
assumption appears to be very restrictive when junction conditions are considered. At least for the modeling 
performed, we have found that combinations of significant physical variables (equations (32) and (35) among 
others) are forced to maintain the first order approximation and, consequently, the most plausible astrophys- 
ical scenario we have obtained within this limit (and EoS), are relativistic rotating matter configurations 
surrounded by an extended “atmospheres”. This is more evident from Tolman Vl-like model (plates C-1, 
C2 and C-3 in figure 7) where, in spite of its singularity, it could represent an object having hydrodynamic 
densities p ~ (10 — 20)/Oo at a core 0 < r ^ 10 Kms. with a thinner matter distribution {p ^ 10^® gr/cm^) 
prolonged to the outer mantle 10 Kms. ^ r ^ 40 Kms. Schwarzschild-like and Tolman IV-like models 
(plates A-1, through B-3 in the same figures) provide the same range for hydrodynamic densities, but for 
configurations having radius 4 times greater than typical ones for neutron stars. We have also found that 
all the obtained models experiment differential rotation, i.e. the core rotates much faster than the envelope. 
But more important than this, is the effect of the dissipation on the orbital velocity. From figures 8, 9, 10 
and 11 it can be appreciated that the more diffusive the model is, the slower it rotates. This effect seems to 
be independent of closure relation and the EoS considered. Notice that we have mainly displayed here figures 
related to the tangential velocity ujz ■ For the complete set of figures for all physical variables corresponding 
every EoS in each scenario and details of some of the calculation, the interested reader is referred to the 
website 

http;//webdelprofesor.ula.ve/ciencias/nunez/EddintonFactor/OtherFiguresVEddingtonFac.html. 
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9 Appendix 

9.1 Einstein equations for slowly rotating matter configuration 

Denoting differentiation with respect to m, r and 9 are denoted by subscripts 0,1 and 2, respectively, Einstein 
equations for slowly rotating matter configuration can be written as: 

* - Guu • 


1-— I ^^2 _ I [p + Pfl ~ +{2u}x + V}T + uj1{P-\-V)\ 


„ / 2 m \ ~ m 22 ni 2 „ „ ~ . 2 m \ ^ 6 m 2 ^ 

2 1-mi--^ cot 6» - 2e-^>^mo + 31-/32- -/32 , 

\r r r \ r J r 


(69) 


• SttTiij- — Gur • 


I + UJa 


[p + PR- - ujj;{P+ V - T)] = 2mi - P 22 - 132 cot (6») - /3| 


(70) 


® 8 V1 -/ 'pj' • 


2..(l-f) (^)[p + P + p. + 7^-2^]=/3i, 


(71) 


• SirTeg = Gee : 


47rr^ [2P — 7^ + p/{] = 2/32 cot 0 — mur — 2e ^^r^/3oi — 6/3imir 


(72) 


2m 


+3/3ir + /3| + 1-(4/3ir + 2/3iir - /3i) r , 


• 87rr„e = G„e : 


2m\ 


0 = I 1-^ I (r/321 + 4r/3i/32 - P 2 ) - ^21 - 2/32mi + — - e ^^r/3o2 - 4/3im2 , (73) 

\ r y r 


• STrTre = Gre : 


0 = /321 -/32 , 

r 


(74) 


• STrTe^ = Ge 


( 2?Ti ^ r~j 

1-^ (4/32/3i - /32i) + m 2 (1 - 4/3ir) - rm 2 i 


2/32 (mir - m^ - e (/32i - /3 o 2) - 2/326 (/3i - /3o) , 


(75) 
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S7rTj-(i) — Grtp • 
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1 - —) {2J^ + V-3pr)V 


~ / 2m 
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2m 


{^x + 1) 


-2c^, ( F + ^ - 1 


(76) 


{iV-pR - 2T) ( 0 ., + 1 - v^l-^x)] } 

= a 2r (/3i + /3o) - 1 + (/3ii - /3oi) + (^1 - 2mi^ + e^^/32 (/32 - 3 cot G) + e^^/32 


• S-TTruai = G 


U<^ — ^ucf) • 


87rr 


2sin0-\/l — a;2 


1-— J (2wa;J^ — 7^ + 3pfl) I?+ 

/ 2m 


[2 {p + Pr — P) — LOx {P PV — P)] + 


Wj; + 1 

+ 'P ~ .7^)1 + „ • n 2^ 

2a;a;Sin0 (1 — uj^) L 

+ x/r^[2.F + u;, (W-pr)]]} 

2m 


iui — 1 


1 - 


P Pr — ^ + 1) ^ 


2ci;a; [p P pft P P — 2P) + 2P (^LOx H“ 1) 


a<ll- 


P (4/3o/3i - 4/3^ + /3oi) - (Po - 3/3i) r - l] 


(77) 


(/3o - 3/3i) (^1 - 2mi^ - moi + mn - ^ (4/317- - 3) + ^1 - (^1 - 2mi^ 


+ e (/3oi - /3oo) - -^ j (3/3| + 6/32 cot 0 + e ^^/3iir^) 

g2/3 ^ ^ 'I 

H—— (^m 22 + 6/32 m 2 + 3m2 j - /32 (/32 - 3 cot 9) - P 22 ? ■ 


9.2 Einstein equations in terms of the flux and Eddington factors 

Einstein equations for slowly rotating matter configuration can be written in terms of the flux and Eddington 
factor as: 
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9.3 The third surface equation 

As we have stated in Section 6, the third surface equation for the slowly rotating Tolman-IV-like model can 
be written as 

q = 0 i + + r 

where the expression for the coefficients 0, $, and F in terms of the Surface Variables are: 

0 = {2n'^F {un^F + 22F + n-30) -n (149P^ + 146F - 37tP + 3 - 7t) + 92F^ - 88F - 4} 
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$ = - {2n^F‘^{mF - 40 + 22F + Tt - 47) + 3 + Tt 

2 FQ 

+0(440F - 149^3 - + 257^^ - 106F - 2) + F(92F - 157F - + 62)} 

with 

U =0F [40F (O + 4) - 220 - 63F + 62] + O + 2F (23F - 22) - 2, 

Q =0F [40F (OF - O + 4F) - (38F - 22) O - (63F - 125) F - 61] - O + F {AQF'^ - 90F + 42) + 2 
and 

n =\/l + 22F - 20OF - 23F2 + 20OF2 + 402F2; 
finally 

r = ^ cfeO^'^) 

where 

Co = 967 rA 3 F( 2 F - 1)(F - 1) 

Cl =2AF{{-8FTL + 6F^ + ATL - 9 + 3F)n + 108F - 267F^ + ISOF^ + 9 - ATL{iF - 2F^ - 1 ) 

+ 7:A^{-576F^ + 864F - 288)) 

C 2 =( 6 F^ + 24ylF^ + 60AF + 12 F - 3 - 15F^ - 8AAF^)n + ttA^F{2112F'^ - 3168F + 1056) - 588y4F^ - 48AF 
- 888 AF^ + 152AAF^ - 3 - 16TLF‘^A{2F - 1) + F{66F^ + 117F - 21 - 159F2) 

C 3 = 6 (F - 1)(-16.4F2 + 2 F + 8AF - l)n + 6 (F - 1){80AF^ - 3AF^ + A9F'^ - ASAF"^ - 14F 
- 96nA^F{2F + 1) - 4AF - 1) 

C 4 = -12F(2F - 1)(F - l)n + 12F(2F - 1)(F - 1)(5F + 8 ^F - 6 ) and C 5 = 24F2(+1 + 2F^ - 3F) 
with 

r =^-0 + l and F = 4^712(20-l)e. 
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Figure 1: The vs / for the different closure relations. L-E^ Lorentz-Eddington; Maximum 

Packing; Mi—s-Mmer^o; W—^Bowers-Wilson; 'M.-C—^Janka (Monte Carlo); Ij-P^Levermore-Pomraning. 
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Figure 2: Variable Flux factor f = f (^x = .^yj^gre /i 

which is considered the Schwarzschild-like and Tolman IV-like models, and changing /2 
for the Tolman Vl-like models. 


fcore = 0.902 
/core = 0.952 
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Figure 3: Profiles of orbital velocity ujz x 10“®c, for the Schwarzschild-like (plates A-1 thought A-3), Tolman 
IV-like (plates B-1 thought B-3) and Tolman VTlike (plates C-1 thought C-3) at three distinct times 
u = 10, 30, 50. The profiles in each plate correspond to a variable Flux factor 

fLEsurface + fcore) 


IlE = fhE 




(l + 
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A-1 A-2 A-3 





r/ao r/ao r/aQ 


Figure 4: Profiles of orbital velocity ujz x 10“®c, corresponding Minerbo closure relation. They are represented 
in A-1 thought A-3 (Schwarzschild-like); B-1 thought B-3 (Tolman IV-like) and plates C-1 thought C-3 
(Tolman Vl-like) at three distinct times u = 10, 30, 50. The profiles in each plate correspond to a variable 

(l -I- 




Flux factor fMi = /mz (x = 
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A-1 A-2 A-3 
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Figure 5: Profiles of orbital velocity ujz x 10“®c, corresponding Janka {Monte Carlo) closure relation. They 
are represented in A-1 thought A-3 Schwarzschild-like); B-1 thought B-3 (Tolman IV-like) and plates C-1 
thought C-3 (Tolman VTlike) at three distinct times u = 10,30, 50. The profiles in each plate correspond to 

. , , r 1 p r ( T \ (fMCsurface T C /core) 

a variable Flux factor Jmc = Jmc J = --' 
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A-1 
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r/aQ 


C-2 



A-3 



B-3 



C-3 



Figure 6: Profiles of orbital velocity lJz x 10“®c, corresponding Maximum Packing closure relation. They 
are represented in A-1 thought A-3 (Schwarzschild-like); B-1 thought B-3 (Tolman IV-like) and plates C-1 
thought C-3 (Tolman VTlike) at three distinct times u = 10,30, 50. The profiles in each plate correspond to 


a variable Flux factor Jmp 



{fMPsurface + ^ /core) 

(l -I- 
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Figure 7: Profiles of hydrodynamic density, p x 10^^ gr/cm^, corresponding to Schwarzschild-like and 
for Tolman IV-like are represented in plates A-l thought A-3 and B-1 thought B-3, respectively. Tol- 
man Vl-like models are displayed in plates C-1 thought C-3 as p x 10^® gr/cm^. The various flux 
factors are Jle = 0.426,0.550,0.750,0.850 for Schwarzschild-like and the Tolman IV-like models and 
Ile = 0.850, 0.900, 0.930 for the Tolman VTlike. The retarded timesdisplayed are u = 10,30, 50. 
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Figure 8: Profiles of orbital velocity ujz x 10“®c, corresponding to Schwarzschild-like for Tolman IV-like and 
Tolman Vl-like are respresented in plates (A-l thought A-3), (B-1 thought B-3) and (C-1 thought C-3), 
respectively. The various constant flux factors are Jle = 0.426,0.550,0.750,0.850 for Schwarzschild-like 
and the Tolman IV-like models and = 0.850,0.900,0.930 for the Tolman Vl-like. The retarded times 
displayed are u = 10,30, 50. 
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Figure 9: Profiles of orbital velocity corresponding to Minerbo closure relation. They are represented in A-1 
thought A-3 (Schwarzschild-like); B-1 thought B-3 (Tolman IV-like) and plates C-1 thought C-3 (Tolman 
Vl-like) with x 10“^c. The constant flux factors are /mi = 0.428,0.550,0.750, 0.850 for Schwarzschild-like 
and the Tolman IV-like models and /mi = 0.910,0.930,0.940 for the Tolman Vl-like with x 10“^c. The 
retarded times displayed are u = 10,30, 50. 
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Figure 10: Profiles of orbital velocity corresponding to Janka {Monte Carlo) closure relation. They are 
represented in A-1 thought A-3 (Schwarzschild-like); B-1 thought B-3 (Tolman IV-like) and plates C-1 
thought C-3 (Tolman Vl-like) with x 10“^c. The constant flux factors are /mc = 0.428,0.550,0.750,0.850 
for Schwarzschild-like and the Tolman IV-like models and /mc = 0.910, 0.930, 0.940 for the Tolman Vl-like 
with ujz X 10“^c. The retarded times displayed are u = 10, 30, 50. 
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Figure 11: Profiles of orbital velocity corresponding to Maximum Packing closure relation. They are repre¬ 
sented in A-1 thought A-3 (Schwarzschild-like); B-1 thought B-3 (Tolman IV-like) and plates C-1 thought 
C-3 (Tolman Vl-like) with xl0“^c. The constant flux factors are /mp = 0.428,0.550,0.750,0.850 for 
Schwarzschild-like and the Tolman IV-like models and /mp = 0.910,0.930,0.940 for the Tolman Vl-like 
with uj^ xl0“^c. The retarded times displayed are u = 10, 30, 50. 
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